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TECHNICAL NOTE 1«)46 


A COMPARATIVE ANALYSIS OF THE PERFORMANCE OF 

LONG-RANGE HYPERVELOCITY VEHICLES 1 

By Alfred J. Eggers, Jr., H. Julian Allen, 
and Stanford E. Neice 


SUMMARY 


Long-range hypervelocity vehicles are studied in terms of their 
motion in powered flight, and their motion and aerodynamic heating in 
unpowered flight. Powered flight is analyzed for an idealized propulsion 
system which rather closely approaches present-day rocket motors. Unpow- 
ered flight is characterized by a return to earth along a ballistic, skip, 
or glide trajectory. Only those trajectories are treated which yield the 
maximum range for a given velocity at the end of powered flight. Aero- 
dynamic heating is treated in a manner similar to that employed previously 
by the senior authors in studying ballistic missiles (NACA TN 4o47 ) > 
with the exception that radiant as well as convective heat transfer is 
considered in connection with glide and skip vehicles. 

The ballistic vehicle is found to be the least efficient of the 
several types studied in the sense that it generally requires the highest 
velocity at the end of powered flight in order to attain a given range. 
This disadvantage may be offset, however, by reducing convective heat 
transfer to the re-entry body through the artifice of increasing pressure 
drag in relation to friction drag - that is, by using a blunt body. Thus 
the kinetic energy required by the vehicle at the end of powered flight 
may be reduced by minimizing the mass of coolant material involved. 

The glide vehicle developing lift-drag ratios in the neighborhood 
of and greater than 4 is far superior to the ballistic vehicle in ability 
to convert velocity into range. It has the disadvantage of having far 
more heat convected to it; however, it has the compensating advantage 
that this heat can in the main be radiated back to the atmosphere. Con- 
sequently, the mass of coolant material may be kept relatively low. 

The skip vehicle developing lift-drag ratios from about 1 to 4 is 
found to be superior to comparable ballistic and glide vehicles in con- 
verting velocity into range. At lift-drag ratios below 1 it is found to 
be about equal to comparable ballistic vehicles while at lift-drag ratios 
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above 4 it is about equal to comparable glide vehicles. The skip vehicle 
experiences extremely large loads, however, and it encounters most severe 
aerodynamic heating. 

As a final performance consideration, it is shown that on the basis 
of equal ratios of mass at take-off to mass at the end of powered flight, 
the hypervelocity vehicle compares favorably with the supersonic airplane 
for ranges in the neighborhood of and greater than one half the circum- 
ference of the earth. In the light of this and previous findings, it is 
concluded that the ballistic and glide vehicles have, in addition to the 
advantages usually ascribed to great speed, the attractive possibility of 
providing relatively efficient long-range flight. 

Design aspects of the glide vehicle are touched on briefly. It is 
argued from considerations of motion and heating that vehicles of this 
type which fly at hypersonic speeds to impact with the earth's surface 
might profitably consist of blunt-nosed bodies of revolution stabilized 
by a conical flare at the base and controlled by deflectable sections of 
the afterbody. In the event that wings are necessary to provide accept- 
able low-speed characteristics, it is indicated that they should have 
highly swept, rounded leading edges in order to alleviate the local 
heating problem with minimum drag penalty. 


INTRODUCTION 


It is generally recognized that hypervelocity vehicles are especially 
suited for military application because of the great difficulty of defend- 
ing against them. It is also possible that for long-range operation, 
hypervelocity vehicles may not be overly extravagant in cost. A satellite 
vehicle, for example, can attain arbitrarily long range with a finite 
speed and hence finite energy input. E. Sanger was among the first to 
recognize this favorable connection between speed and range (ref. l) and 
was, with Bredt, perhaps the first to exploit the speed factor in design- 
ing a long-range bomber (ref. 2) . This design envisioned a rocket-boost 
vehicle attaining hypervelocity speeds at burnout and returning to earth 
along a combined skip-glide trajectory. Considerable attention was given 
to the propulsion and motion analysis; however, little attention was given 
to what is now considered to be a principal problem associated with any 
type of hypersonic aircraft, namely that of aerodynamic heating. In 
addition, the category of expendable vehicles, perhaps best characterized 
by the ballistic missile, was not treated. 

Since the work of Sanger and Bredt there have been, of course, many 
treatments of long-range hypervelocity vehicles in which the propulsion, 
motion, and heating problems have been studied in considerable detail. 
However, these analyses have been devoted in the main to particular designs 
and are not intended to reveal, for example, the relative advantages and 
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disadvantages of ballistic-, skip-, and glide-type vehicles. Furthermore, 
it appears that the extent to which these vehicles can compete on a simple 
efficiency basis with lower speed aircraft of either the expendable or 
nonexpendable type has not been well established. 

It has therefore been undertaken in the present report to make a 
comparative analysis of the performance of hypervelocity vehicles having 
ballistic, skip, and glide trajectories. An idealized propulsion system, 
whose performance is rather closely approached by present-day rocket 
motors, is assumed. The motion analysis is simplified by treating, for 
the most part, only optimum trajectories yielding the maximum range for 
a given initial kinetic energy per unit mass in the unpowered portion of 
flight. Aerodynamic heating is treated in a manner analogous to that 
employed by the senior authors in studying ballistic missiles (ref. 3) 
with the exception that radiant heat transfer, as well as convective heat 
transfer, is considered in the treatment of glide and skip vehicles. The 
efficiencies of these vehicles are compared with supersonic aircraft with 
typical air-breathing power plants. 


ANALYSIS 

General Considerations 


In the following analysis of long-range hypervelocity vehicles, only 
flight in planes containing the great circle arc between take-off and 
landing is considered. The flight is thought of in two phases: (a) the 

powered phase in which sufficient kinetic energy, as well as control, is 
imparted to the vehicle to bring it to a prescribed velocity, orientation, 
and position in space; and (b) the unpowered phase, in which the vehicle 
travels to its destination under the influence of gravity and aerodynamic 
forces. 

The analyses of motion and aerodynamic heating during unpowered flight 
will, of necessity, differ widely for the several types of vehicles under 
consideration. On the other hand, motion in the powered phase is con- 
veniently treated by a method common to all vehicles. The study of powered 
flight and its relation to range is therefore taken as a starting point 
in the analysis. The mathematical symbols employed in the analysis are 
listed in Appendix A. 


Powered Flight and the Breguet Range Equation 


In this part of the study, the following simplifying assumptions are 
made: (a) aerodynamic heating can be neglected on the premise that high 
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flight speeds are not attained until the vehicle is in the rarefied upper 
atmosphere; 2 (b) sufficient stability and control is available to provide 
proper orientation and positioning of the vehicle in Space; (c) the dis- 
tance traveled while under power is negligible by comparison to the over- 
all range; and finally, (d) the thrust is very large compared to the 
retarding aerodynamic and gravity forces. In terms of present-day power 
plants, the last assumption is tantamount to assuming a rocket drive for 
the vehicle . 

The velocity at burnout of the first stage of a multistage rocket 
(or the final velocity of a single-stage rocket) can then be expressed 
as (see, e.g., ref. 4): 


gl 

- Vi 



( 1 ) 


where the initial velocity is taken as zero. In this expression, m^ and 
mf represent the mass of the vehicle at the beginn ing and ending of 
first-stage flight, and = Vf^/Vg where Vg = V gr Q = 25,930 feet 

per second is the satellite velocity at the surface of the earth. The 
coefficient g is the acceleration due to gravity and is, along with 
the specific impulse I, considered constant in this phase of the anal- 
ysis. The final velocity of the vehicle at the end of the N stages of 
powered flight can be -expressed as 


Vf = V N = In 




( 2 ) 


where the initial mass of any given stage differs from the final mass of 
the previous stage by the amount of structure, etc., jettisoned. 

Now let us define an equivalent single-stage rocket having the same 
initial and final mass as the N-stage rocket and the same initial and final 
velocity. There is, then, an effective specific impulse defined by 


2 This assumption is in the main permissible. A possible exception 
occurs, however, with the glide vehicle for which heat-transfer rates 
near the end of powered flight can be comparable to those experienced in 
unpowered gliding flight. 
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whereby equation (2) can be written as 


v f - ln (s?) C*) 

The effective specific impulse I e is always somewhat less than the 
actual specific impulse , but for an efficient design they are very nearly 
equal. Throughout the remainder of the analysis the effective impulse 
I e will be used. 

Equation (4) might be termed the "ideal power plant" equation for 
accelerated flight because, when considered in combination with the 
assumptions underlying its development, attention is naturally focused 
on the. salient factors leading to maximum increase in velocity for given 
expenditure of propellant. Thus the thrust acts only in overcoming inertia 
forces, and the increase in vehicle velocity is directly proportional to 
the exhaust velocity (gl) for the propellant. 

Now we recognize that an essential feature of the hypervelocity 
vehicles under study here is that they use their velocity (or kinetic 
energy per unit mass) to obtain range. For this reason, equation (4) also 
constitutes a basic performance equation for these vehicles because it 
provides a connecting link between 'range requirements and power-plant 
requirements. 

In addition to comparing various types of hypervelocity vehicles, 
our attention will also be focused upon comparison of .these vehicles with 
lower speed, more conventional types of aircraft. For this purpose it is 
useful to develop an alternate form of equation (4). We observe that the 
kinetic energy imparted to the vehicle is 



2 


This energy is equated to an effective work done, defined as 
of the range traveled and a constant retarding force. (Note 
useful kinetic energy at the end of powered flight is zero.) 
is termed the "effective drag" D e . Thus 


the product 
that the 
This force 
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DeR 


= | m f v f 2 


(5) 


where R is flight range measured along the surface of the earth. 
Similarly, we may define an "effective lift” L e , equal to the final 
weight of the vehicle 

L e = W f = m f g 


from which it 


follows that equation (5) may 


R.ii 


V* 


P/e 2g 


Toe written as 


(6) 


where (L/D) e is termed the "effective lift-drag ratio." Combining 
equations ( 4) and ( 6) , we obtain 


R 



I I e V e ln 
e 



(7) 


where 



( 8 ) 


and represents an "effective" flight velocity of the vehicle. Equa- 
tion (7) will prove useful in comparing hypersonic vehicles with conven- 
tional aircraft because of its analogy to the Breguet range equation, 

E - D W ln (®?) <9) 

It will also prove useful to have equation (7) in the dimensionless form 
obtained by dividing through with r 0 , the radius of the earth. In this 
case we have 



where $ is the range in radians of arc traversed along the surface of the 
the earth. 


Motion in Unpowered Flight 


Ballistic trajectory .- In studying the motion of long-range vehicles 
in this trajectory, advantage is taken of the fact that the traverse 
through the earth's atmosphere generally forms only a small part of the 
total trajectory. Therefore, the deflection and deceleration encountered 
in the re-entry phase (discussed in detail in ref. 3) are neglected in 
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the computation of the total range and rotation of the earth is neglected 
in this and al 1 other phases of the analysis. With the added simplifica- 
tion that the contribution to range of the powered phase of flight is 
negligible, the ballistic trajectory becomes one of Kepler's planetary 
ellipses, the major axis of which bisects the total angle of arc $ 
traveled around the earth (see sketch). For the trajectories of interest 
here (V f <l), the far focus of the ellipse is at the mass center of the 
earth. For purposes of range computation, then, the ballistic vehicle 
leaves and returns to the earth's surface at the same absolute magnitude 
of velocity and incidence (see sketch). 



where the angle of incidence Of is considered positive. In order to 
determine the optimum trajectory giving maximum range for a given velocity 
Vf, equation (ll) is differentiated with respect to Of and equated to 
0, yielding 

Vf = — = 1 - tan 2 0f 

V 

$ = ^P- = it - hef 



( 12 ) 
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Equations (ll) and (12) have been employed to determine velocity as a 
function of incidence for various values of range and the results are 
presented in figure 1. The "minimum velocity line" of figure 1 cor- 
responds to the optimum trajectories (eqs. (12)). 

The effective lift-drag ratios can easily be calculated for optimum 
ballistic vehicles using equation (6) in combination with the information 
of figure 1. The required values of (L/d) e as a function of range are 
presented in figure 2. 

Skip trajectory .- This trajectory can be thought of as a succession 
of ballistic trajectories, each connected to the next by a "skipping 
phase" during which the vehicle enters the atmosphere, negotiates a turn, 
and is then ejected from the atmosphere. The motion analysis for the 
ballistic missile can, of course, be applied to the ballistic phases of 
the skip trajectory. It remains, then, to analyze the skipping phases 
and to combine this analysis with the ballistic analysis to determine 
over-all range. 

To this end, consider a vehicle in the process of executing a skip 
from the atmosphere ( see sketch) . 



W 


Earth's surface 
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The parametric equations of motion in directions perpendicular and 
parallel to the flight path s are, respectively. 


Cl A - mg cos 0 = 5^- 
L 2 r c 

-Cn A + mg sin 0 = m — 
V 2 dt 


( 13 ) 


where r c is the local radius of curvature of the flight path, 0 is the 
local inclination to the horizontal (positive downward) , p is the local 
air density, and Cl and Cp are the lift and drag coefficients, respec- 
tively, based on the reference area. A, .of the aircraft. 


In the turning process, aerodynamic lift must obviously predominate 
over the gravity component, mg cos 0. By analogy to the atmospheric 
re-entry of ballistic missiles (see ref. 3)j aerodynamic drag generally 
predominates over the gravity component, mg sin 0. Moreover, the inte- 
grated contribution to velocity of this gravity component during descent 
in a skip is largely balanced by an opposite contribution during ascent. 

For these reasons we will idealize the analysis by neglecting gravity 
entirely. This approach is analogous to the classical treatment of impact* 
problems in which all forces exclusive of impact forces (aerodynamic forces 
in this case) are neglected as being of secondary importance. Gravity is 
shown to be of secondary importance in figure 3 where the trajectory 
results obtainable from equations ( 13 ) and (l^) are presented for the first 
skipping phase of an L/D = 2, $ = 1 skip missile. 

With gravity terms neglected, equations ( 13 ) reduce to 


1 C L pV 2 A = -mV 2 — 

2 ds 


- 1 <*»*** - “ § 


( 14 ) 


where de/ds = — — to the accuracy of this analysis. 


Now we assume an isothermal atmosphere, in which case 


p = p 0 e'Py (15) 

where p 0 and |3 are constants, and y = (r - r 0 ) is the altitude from 
sea level (see ref. 3 for discussion of accuracy of this assumption). 
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Noting that dy/ds = -sin 0, we combine the first of equations (l4) with 
equation (15) to yield 


ClPq a 

2m 



sin 0 d0 


(16) 


This expression can be integrated to give 


e“Py = cos 0 - cos 0 e n (17) 

23m 

where p is taken as zero at the altitude corresponding to the effective 
"outer reach" of the atmosphere. Equation (17) points out an important 
feature of the skip path; namely, cos 0 is a single -valued function of 
altitude. Since 0 proceeds from positive to negative values, it is 
evident that 

8en n-i =-8ex n 


where the subscripts en and ex refer to atmospheric entrance and exit 
conditions, respectively, and the numbers n - 1 and n refer to suc- 
cessive ballistic phases of the trajectory. Now since 


dV _ v dY = 1 dV 2 
dt ds 2 ds 


equations (l4) may be combined to obtain 

1 d 7 f = _Vf_ M (19) 

2 ds L/D ds 

which, for constant L/D, can be integrated to yield 

8 ex n - e en n . 1 ( 2 0) 

v ex n _ „ L/D 

T T — “ 


With the aid of equation (l8 ) , 


this expression may he written 
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v ex n 

^ en n-i 


20en n-i 


= e 


L/D 


( 21 ) 


which relates the velocities at the beginning and end of a skip to the 
lift-drag ratio and the entrance angle of the vehicle to the earth's 
atmosphere. From equation (18) it follows further that the entrance 
angle for each skip in the trajectory is the same, so that 


0 en n - 0 en n _ 1 

and hence equation (2l) becomes 


Vex n 

V en n -i 


- 


29f 

L/D 




We now combine this result of the skip analysis with that of the 
ballistic analysis to obtain the total flight range. From equation (ll), 
the range of the nth ballistic segment of the trajectory is 


cp n = 2 tan" 1 


sin 0j*cos df 

(wX - cos2ef 


(23) 


Consistent with the idealization of the skipping process as an impact 
problem, we neglect the contribution to range of each skipping phase so 
that the total range is simply the sum of the ballistic contributions. 
From equations (22) and (23) this range is then 


0 



-z 

n=n 


oo 



n=o 



(24) 
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From this expression we see that for any given velocity at the end 
of powered flight there is a definite skipping angle which maximizes the 
range of an aircraft developing a particular lift-drag ratio. These 
skipping angles have heen_obtained with the aid of an IBM CPC , and the 
corresponding values of V f as a function of range for various l/d are 
presented in figure 4. Corresponding values of (L/D) e have been obtained 
using equation (6) and the results are shown in figure 5» 

Glide trajectory.- The trajectory of the glide vehicle is illustrated 
in the accompanying sketch. As in the previous analyses, the distance 
covered in the powered phase will be neglected in the determination of 
total range . 



The parametric equations of motion normal and parallel to the direc- 
tion of motion are the relations of equations ( 13 ) rewritten in the form 


L - mg cos 0 


mV 2 


> 


-D + mg sin 0 


m 


dV 

dt 


(25) 


Under the assumption of small inclination angle 0 to the horizontal 
(thus cos 0 = 1, sin 0=0), constant gravity acceleration (i.e., — ~l), 
and noting the following relations 
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dV = dV = 1 dvf 

dt ds 2 ds 

_1_ = d( ^ - e) 
r c ds 

d]jf _ cos 9 „ _1_ 
ds r r Q 


( 26 ) 


equations ( 25 ) can "be written in the forms 


t 2 d£ _i_ mV 2 

L = mV — + mg 

ds & r Q 


(27) 


D = - — m + mg 9 
2 ds 

Dividing the first of equations (27) "by the second yields the following 
differential equation 


g 



\ fl L dV 2 
/ + \2 D ds 


+ V c 




( 28 ) 


But, as is demonstrated in Appendix B , the terms ^ gO and V 2 may- 
be neglected so that equation ( 28 ) reduces to 


dV 

ds 


,(L/D) 


V 2 + 


2g 

l/d 


= 0 


(29) 


Since 

, r 2 

V S = ^o 

equation ( 29 ) can he integrated for constant ^ to give the velocity 
in non ’.imensional form as 
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This expression gives velocity as a function of range for what Sanger 
(ref. 2 ) has termed the equilibrium trajectory - that is, the trajectory 
for which the gravity force is essentially balanced by the aerodynamic 
lift and centrifugal force, or 


L 

W 



(31) 


It follows from equation (25) that velocity can be expressed in the form 


V 2 


1 

ClAV S 2 P 

2 mg 


(32) 


Now it is intuitively obvious that as the maximum range is approached, 
L/w —>1 and hence V 2 becomes small compared to one (see eq. (3l))* 

In this event it follows from equation (30) that the maximum range for 
the glide vehicle is given by 


$ 


_R_ 

r o 


1 

2 





(33) 


The relation between velocity and range has been determined with 
equation ( 33 ) for various values of l/d and the results are presented 
figure &• Corresponding values of (l/d )0 have been obtained using 
equation ( 6 ) and are presented in figure 7 . 

These considerations complete the motion analysis and attention is 
now turned to the aerodynamic heating of the several types of vehicles 
under consideration. 


Heating in Unpowered Flight 


General considerations.- Three aspects of the aerodynamic heating 
of hypervelocity vehicles will be treated here; namely, 

1. The total heat input 

2. The maximum time rate of average heat input per unit area 

3. The maximum time rate of local heat input per unit area 
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Total heat input is, of course, an important factor in determining 
over-all coolant weight, whether the coolant he solid (e.g., the struc- 
ture), liquid, or gas, or a combination thereof • The maximum time rate 
of average heat input per unit area can determine peak average flow 
rates in the case of fluid coolants and may dictate over-all structural 
strength in the event that thermal stresses predominate. 

Excessive local heating is, of course, a serious problem with hyper- 
velocity vehicles. This problem may vary depending upon the type of the 
vehicle. Thus, for the ballistic vehicle, an important local "hot spot" 
is the stagnation region of the nose, while for the skip or glide vehicle 
attention may also be focused on the leading edges of planar surfaces 
used for developing lift and obtaining stable and controlled flight. In 
this analysis attention is, for the purpose of simplicity, restricted to 
the "hot spot" at the nose. In particular, we consider the maximum time 
rate of local heat input per unit area because of its bearing on local 
coolant flow rates and local structural strength. 

It is undertaken to treat only convective heat transfer at this 
stage of the study. As will be demonstrated, radiant heat transfer from 
the surface should not appreciably influence convective heat transfer 
to a vehicle. Therefore, alleviating effects of radiation are reserved 
for attention in the discussion of particular vehicles later in the paper. 
This analysis is further simplified by making the assumptions that 

1. Effects of gaseous imperfections may be neglected 

2. Shock-wave boundary-layer interaction may be neglected 

3. Prandtl number is unity 

k. Reynolds analogy is applicable 


These assumptions are obviously not permissible for an accurate quanti- 
tative study of a specific vehicle. Nevertheless they should not inval- 
idate this comparative analysis which is only intended to yield informa- 
tion of a general nature regarding the relative merits and problems of 
different types of vehicle (see ref. 3 for a more complete discussion 
of these assumptions in connection with ballistic vehicles). 


In calculating convective heat transfer to hypervelocity vehicles, 
the theoretical approach taken in reference 3 for ballistic vehicles is, 
up to a point, quite general and can be employed here. Thus, on the 
basis -of the foregoing assumptions, it follows that for large Mach numbers, 
the difference between the local recovery temperature and wall tempera- 
ture can be expressed as 


l 



(3*0 


(T r - T w ) 
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where Tr is the recovery temperature, T w is the wall temperature, Cp 
is the specific heat at constant pressure, and the subscript l denotes 
local conditions at a point on the surface. 

It is clear, however, that the walls of a vehicle should he maintained 
sufficiently cool to insure structural integrity. It follows in this case 
that, the recovery temperature at hypervelocities will be large by compar- 
ison to the wall temperature and equation (34) may be simplified to read 



(35) 


To the accuracy of this analysis, then, the convective heat transfer is 
independent of wall temperature. Therefore, as previously asserted, 
radiant heat transfer should not appreciably influence convective heat 
transfer and the one can be studied independently of the other. 

Now, according to Reynolds analogy, the local heat-transfer coef- 
ficient h i is, for a Prandtl number of unity, given by the expression 


h l ~ 2 C Fi C Vi P l V l 


(36) 


where Cp is the local skin-friction coefficient based on conditions 

just outside the boundary layer. With the aid of equations (35) and ( 36 ) 
the time rate of local heat transfer per unit area, 


f| = h z (T r - T w ) 1 


(37) 


can be written as 


dH _J1 
dt 4C p 


cf 2 c Pz p z v z 


(38) 


Equation ( 38 ) can be integrated over the surface of a body to yield the 
time rate of total heat input as follows 


dQ 

dt 


dS = 7 pV 3 C F ’S 
dt 4 


S 


(39) 
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vherein 



is set equal to 


C P 


and 


C\F f - 



p 7 V Z 

— — dS 
P V 


m 


The parameter Cp* is termed the "equivalent skin-friction coefficient" 
and will be assumed constant for a particular vehicle. From equation ( 39 ) 
we can obtain two alternate forms which will prove useful; namely, the 
altitude rate of total heat input defined by (note that dy is negative 
for dt positive) 


IQ 1 dQ _ pV 2 Cp t S 

dy V sin Qf dt 4 sin 0 f 


and the range rate of total heat input defined as 


dQ 1 dQ = P y2 CF tS 

d(r Q q>) V cos ®f dt h cos 0 f 


(42) 


The total heat input may be obtained by integration of equations (39), 
(4l) or (42), depending upon the particular variable used. 

The time rate of average heat input per unit area may be obtained 
from equation ( 39 ) as 


dHav 

dt 


= 1 dQ = 1 
S dt 4 


pV°Cx 


(43) 


Consider next the local convective heat transfer in the region of 
the nose. The time rate of local heat input per unit area was determined 
in reference 3 under the assumptions that viscosity coefficient varies as 
the square root of the absolute temperature, and that flow between the 
bow shock wave and the stagnation point is incompressible. In this case 
it was found that 
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^ = 6.8x10-®/! V 3 (44) 

where a is the radius of curvature of the body at the stagnation point. 

With these relations we are now in a position to study the heating 
of the several types of vehicles of interest. 


Ballistic vehicle .- The heating for this case has already been 
analyzed in reference 3. Only the results will be given here. 

The ratio of the total heat input to the initial kinetic energy was 
found to be 


Q 


Cp* S 


c d p 0 a 

Pm sin 6f 


1 „ 2 
- mVf 


m 


2 C D A 

For the "relatively light missile," which is of principal interest here, 


CpPpA 

Pm sin Qf 
e < < 1 


(46) 


and equation (45) reduces to 


Q 


i c F «s 


i mV f 2 2 C D A 

o 1 


(47) 


The time rate of average heat input per unit area was found to be 

(48) 


3 c dPqA m _j3y 


dHav 1 dQ. CF'PpVf -py 2Pm sin 0f 

~dt S dt 4 6 


which has the maximum value 


^av\ __ /dH a y\ 

. dt /max V dt /y a 


— mV : p 3 sin 6f 
6e 0j)A 


(*9) 


at the altitude 
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3C d p q a \ 
,23m sin Of) 


(50) 


Equation (49) applies, of course, only if the altitude of occurrence is 
above ground level. If the value of y a is negative then the maximum 
rate will, of course, occur at sea level. 

The time rate of local heat input per unit area to the stagnation 
region of the nose was found to be 


By 3CDPnA -Py 

= 6.8x10 »e' 2 e sin 8 f 
dt. v a f 


(51) 


having a maximum value of 




6.8x10 


r Pm sin df 
3eCporA 



(52) 


occurring at the altitude 


1- 

^b = p 


In 


" 3C d p q a \ 
Pm sin Of) 


(53) 


If the value of y-^ is negative, then the maximum value occurs at 
ground level. 

Skip vehicle .- With the aid of equation (17), the density at any 
point in a given skipping phase is found to be 


P 


Po e 


-Bv 2pm . . 

F,y = (cos 6 - cos 0f) 


(54) 


where it is to be recalled that 0 en = Of. The corresponding velocity 
is obtained by integrating equation (19) for constant L/D, yielding 
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V 


Ven e 


0f-0 

L/D 


(55) 


By substitution of equations (54) and (55) into equation (39), the time 
rate of total heat input at any point in a skipping phase can be expressed 
as follows: 


dQ _ 1 Cf*S Pm 
dt _ 2 CpA l/D 


V en 3 (cos 0 


3(9f-e) 

cos 0f)e 


(56) 


Now, recalling that ds/dt = V, the first of equations (l4) may be com- 
bined with equation (17) to yield 


= -pv( cos 0 - cos 0f) 
dt 


(57) 


Inasmuch as V en = V ex for any ballistic phase, it then follows from 
equation (22) that 


(Ven) D _ (^eJn _ El (58) 

(V en )n = i Vf 

With the aid of equations (57) and (58), equation (56) can be integrated 
to give the total heat input for a given skipping phase. Thus we obtain 


Qn 


1 c F «s 

2 C D A 




(59) 


where n refers to the ballistic phase subsequent to the given skipping 
phase. 


The total heat input for the entire trajectory can be obtained by 
su mm ing up the heat inputs for each separate skipping phase. Performing 
this operation yields 
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= y Qn 

1 mVf 2 L> ^ mVf 2 

2 1 n=l 2 1 


1 C F »S 

2 C D A 




V - (n-1) 

A e 

n=i 


ke t 

L/D 


(60) 


The summation on the right side of equation (60) represents a geometric 
series which can readily be evaluated. The total heat input for the 
entire trajectory then becomes 


Q _ 1 C F* S 
1 mVf 2 2 C D A 


(6i) 


which is identical to the result obtained for the light ballistic missile 

(eq. (47)). 


The time rate of average heat input per unit area is obtained by 
dividing equation ( 5 6 ) with the surface area, thus yielding 


dH av 

dt 


1 C F ' pm 

2 CpA l/D 


3(0f-e) 

V e n 3 ( cos 9 ~ cos 0f)e ^/D 


(62) 


It can be shown that this expression has a peak value at a point in the 
skip, 0 a , given by 


or 



(63) 


From equation (22) it can be concluded that the maximum heat-transfer 
rate will occur in the first skip where V en = Vf; consequently. 
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3(0f-e a ) 
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(64) 


The time rate of local heat input per unit area in the stagnation 
region of the nose is obtained by introducing equations ( 54 ) and ( 55 ) 
into equation ( 44 ) with the following result: 


dH s 

dt 
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(65) 


Equation (65) has a peak value at a point 0 ^ in a skip given by 


(cos 0-^ - cos 0f) 


L/D 

= — sin 0 b 


or 


0 b = tan -1 


6 

L/D 


sin 


cos 0f 



( 66 ) 


It is clear in this case also that the heat-transfer rate will have its 
maximum value in the first skipping phase where the velocities are 
highest. Since V en = V f in the first skip, equation (65) becomes 



6.8x10 


-6 


Pm(L/D) sin 0 b 
3C L Aa 


1/2 


Vf 3 e 


3(0f-0 b ) 

L/D 


(67) 


Glide vehicle .- From equations ( 30 ) and (32), the density at a point 
in the glide trajectory is found to he 
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2mg (1 - Vf 2 ) e J 
C L AV S 2 1 - (1 - Vf 2 ) e J 


( 68 ) 


where 


= _ap_ 

l/d 


By substitution of equations (30) and (68) into equation (39), the t im e 
rate of total heat input can he expressed as 


d£ _ 1 C F >S mgVs , 
dt 2 c d A L/D 



- (1 - v f 2 ) 


-1I/2 


(69) 


Now with the aid of equations (30) and (33), equation (69) can be inte- 
grated over the limits of the glide trajectory to yield the total heat 
input in terms of the initial kinetic energy as 



1 C F *S 

2 "CpA* 


(70) 


which expression is identical with that obtained for the skip vehicle 
(eq. (6l)) and for the light ballistic vehicle (eq. (47)). 

Now the time rate of average heat input per unit area is found by 
dividing equation (69) with the surface area, thus yielding 


flHav _ 1 (V mgVg 
dt 2 C d A l/D 


(1 


Vf 2 


)e L 
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V f 2 ) 


-1I/2 



(71) 


It follows from this expression that the maximum time rate of average 
heat input per unit area is 
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/dH av \ _ _ 1 _ mgVs 
V dt /J a 3 n/3 C D A l / d 


(72) 


at a value J a given by 



(73) 


If J a is taken as a reference value, and equations (71) and (72) expres- 
sed in terms of J a and incremental changes AJ = J - J a , it can easily 
he shown that 


dH a y / dt 
(<ffl av /dt) 

max 



- 2e 


AJ\ 


1/2 


J 


= F a (AJ) 


(74) 


The dependence of F a (AJ) on AJ is shown in figure 8. 

The velocity at which the maximum average heat input rate occurs 
can be obtained by substituting equation (73) into equation (30) yielding 


V 



(75) 


In equation (30) it is seen that the velocity is greatest at the start of 
unpowered flight (i.e., when _J = 0) . Equations (72), (73), and (74) 
apply, therefore, only when Vf > (l//3). 

For cases when Vf S (l/\/3)> the maximum time rate of average heat 
input per unit area will occur at the start of unpowered flight and is 
given by 



'dH av \ 1 C F ' mgV s 

dt J 2 c d a l/d 

J a =o 


(1 - Vf 2 ) Vf 


(76) 
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The maximum time rate of local heat input per unit area in the stag- 
nation region of the nose is found by first substituting equations (30) 
and (68) into equation (44) to obtain 


Es = 6 < 8 x 10 - C ^_ 2 

dt V C L Aa 


vs‘ 


l-(l-V f 2 )e J (1 - Vf 2 )e J 


][< 


ll/2 


(77) 


The maximum time rate is then 


dHA 

'max 


dt /J b 373 


8x10 v s 2 


(78) 


occurring at a value of J b given by 


J b = -in 3(1 - V f 2 ) 


(79) 


With as a reference, it can easily be shown that 

dH s /dt 
( as s /dt) max . 

where 

AJ = J - J b 


AJ 

I * 2 


e AJ ) = F b (AJ) 


( 80 ) 


The dependence of F b (Aj) on AJ is shown in figure 8. 

With reference to equations (30) and (74) it can be seen that the 
maximum time rate of local heat transfer in the stagnation region occurs 
when 



( 81 ) 
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It is apparent then that equations ( 77 ), (78), and (79) apply only when 
Vf >-72/3. For cases where V f <J 2/3, the maximum time rate of local 

heat input per unit area will occur at the start of unpowered flight and 
is given by 


\ dt /max \ dt / j^_ 0 




(82, 


DISCUSSION 

Performance of Hypervelocity Vehicles 


In this study the point of view is taken that the performance of 
long-range hypervelocity vehicles is measured by their efficiency of 
flight. Thus, for example, it is presumed that the advantages (military 
and otherwise) of short time of flight accrue equally to all vehicles. 

The efficiency of flight is perhaps best measured by the cost of 
delivering a given pay load a 'given range - the higher the cost, the 
lower the efficiency. Quite obviously it is far beyond the scope of the 
present paper to actually compute this cost. Rather, then, we adopt a 
more accessible parameter of hypervelocity flight, namely, the initial 
mass of the vehicle, as a measure of cost. In effect, then, the assump- 
tion is made that the higher the initial mass of a vehicle the higher the 
cost and the lower the efficiency. With these thoughts in mind, it is 
constructive to reconsider the basic performance equation (eq. (4)) writ- 
ten in the form 


mf 


mfe 


v f/gl e 


(83) 


This expression clearly demonstrates the roles played by the three factors 
which influence the initial mass of a vehicle required to travel a given 
range. For one thing there is the power plant, and as we would expect, 
increasing the effective specific impulse increases the over-all efficiency 
of flight in the sense that it tends to reduce the initial mass. The 
velocity at burnout influences initial mass by dictating the amount of 
fuel required, and it is not surprising that decreasing the required burn- 
out velocity (e.g., by increasing the L/D of a skip or glide vehicle) 
tends to decrease the initial mass. Finally, we see that the initial 
mass is proportional to the final mass which consists of the pay load, 
structure (and associated equipment), and coolant. If we presume the mass 
the pay load to be some fixed quantity , then the initial mass will 
vary in accordance with this mass of structure and coolant. 
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Wow we assume for comparative purposes that the power plant for one 
vehicle is equally as good as the power plant for another vehicle - that 
is to say I e is a more or less fixed quantity. In this event it is per- 
missible to restrict our attention to two main performance considerations; 
namely, the prescribed motion as it influences the required burnout 
velocity, and the resulting aerodynamic heating as it influences structure 
and coolant. We therefore proceed to discuss the comparative performance 
of long-range hypervelocity vehicles in terms of these considerations. 

Motion . - The dependence of burnout velocity Vf on range was deter- 
mined in the analysis of motion in unpowered flight and the results 
obtained for the several types of hypervelocity vehicle under study were 
presented in figures 1, 4, and 6. Using these results in combination 
with the basic performance equation we have calculated the corresponding 
initial to final mass ratios m±/mf as a function of range. For these 
and subsequent calculations it has been assumed that the rocket power 
plant develops an effective specific impulse of 225 seconds. The results 
of these calculations are presented in figure 9 and we observe that, in 
general, the mass ratios are highest for the ballistic vehicle. The glide 
and skip vehicles have comparable and relatively low mass ratios at lift- 
drag ratios in the neighborhood of 4 and greater. The skip vehicle is 
superior, however, to the glide vehicle at lift-drag ratios in the neigh- 
borhood of 2. From considerations of motion alone, then, we conclude 
that the skip vehicle and the glide vehicle developing lift-drag ratios 
greater than 2 are superior efficiencywise, in the sense of this report, 
to the ballistic vehicle. Let us now determine how these observations 
are modified by considerations of aerodynamic heating. 

Aerodynamic heating .- The analysis has revealed one particularly 
salient factor in regard to the heat transferred by convection to hyper- 
velocity vehicles that expend the majority of their kinetic energy of 
flight in traveling through the earth's atmosphere. This factor is that 
the amount of kinetic energy which appears in the body in the form of 
heat is proportional to- the ratio of friction force to total drag force 
acting on the body (see eqs. (47), (6l), and (70)). With the possible 
exception of the relatively heavy ballistic vehicle (see ref. 3) all of 
the hypervelocity vehicles treated here do expend the major part of their 
kinetic energy in flight. It is, in fact, only by virtue of this expend- 
iture of energy that the skip and glide vehicles achieve long range. From 
the standpoint, then, of reducing the total heat transferred by convec- 
tion, the problem is to determine how the ratio of friction 'force to total 
drag force can be reduced. This matter was discussed in detail In ref- 
erence 3 in connection with ballistic vehicles and it was demonstrated 
that the ratio could be reduced by employing high-pressure -drag (i.e., 
blunt) shapes. It would be most fortunate if this avenue of solution 
were open also to the skip and glide vehicles; however, it is readily 
apparent that such is not the case. This conclusion follows simply from 
the fact that the skip and glide vehicles must develop reasonably high 
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lift-drag ratios to achieve long range. But, as is well known, high 
lift-drag ratios and high pressure drag are incompatible aerodynamic 
properties. Evidently, then, the skip and glide vehicles will be rela- 
tively slender and they will, by comparison to blunt ballistic vehicles, 
be required to absorb large amounts of their kinetic energy of flight m 
the form of heat. On the basis of the calculations of reference 3, it 
does not seem feasible for slender hypervelocity vehicles to absorb and 
retain so much heat (of the order of one-tenth the kinetic energy of 
flight). We are led, therefore, to consider the possibility of radiating 
part or all of this heat back to the atmosphere. 

Let us first consider radiation heat transfer from the surface of a 
glide vehicle. For purposes of simplicity we presume a vehicle conical 
in shape. The base diameter is taken as 3 feet and the weight as 5,000 
pounds. We consider two slender cones which, according to hypersonic 
theory including friction drag, can develop maximum lift -drag ratios o 
4 and 6 (see Appendix C) . We find (see Appendix D) that the L/D - 4 
elide vehicle can radiate heat like a black body at a rate equal to the 
maximum- average convective heat-transfer rate if the surface temperature 
is allowed to rise to about 1500° F. If the vehicle develops a lift- 
drag ratio of 6, then the allowable surface temperature must be increased 
to about 1800° F. These surface temperatures are high; nevertheless they 
are within the range of useful strengths of available alloys. Further- 
more, they can, if necessary, be reduced somewhat by designing a less 
dense vehicle (or, more specifically, a vehicle of lower mg/S, see 
Appendix D) . 

It is indicated, then, that the glide vehicle has the attractive 
possibility of radiating back to the atmosphere a large fraction of the 
heat transferred to it by convection. As a result the mass of coolan 
required to protect the vehicle may be greatly reduced. Just as with 
the ballistic vehicle, however (see ref. 3), it is evident that additional 
means, such as transpiration cooling, may well be necessary to protect 
local hot spots on the surface, like the stagnation region of the nose. 

It is also well to note that the alleviating effects of radiative cooling 
are not limited to the glide vehicle alone, but would apply to any hyper- 
velocity vehicle in level flight. 

We inquire now if the skip vehicle is capable of radiating heat at . 
a rate comparable to the maximum convective heat-transfer rate. For this 
purpose it suffices to confine our attention to the first skip wherein 
the maximum convective heat-transfer rates are encountered (see ej. (64)). 
On the basis of our calculations for glide vehicles developing lift-drag 
ratios of 4 and 6, we conclude that the skip vehicle developing compa- 
rable lift-drag ratios cannot radiate heat at anything like the maximum 
convective rate. This conclusion follows directly from the fact that, 
although the heat absorbed by the skip vehicle in the first skip would 
be about the same as that for the glide vehicle experiencing the same 



NACA TN UO 46 


29 


loss of kinetic energy, the rates of absorption would be far greater for 
the skip vehicle. Hence, the surface temperatures required for radiation 
to offset convection would probably exceed the temperatures' at which 
known materials retain appreciable strength. 

Now the skip vehicle operating at lift-drag ratios in the neighbor- 
hood of 2 will absorb less heat than skip vehicles developing higher 
lift-drag ratios. However, as shown in Appendix D, the former vehicle 
still absorbs more heat than a comparable high-pressure-drag ballistic 
vehicle and it accrues no appreciable advantage by radiation. From the 
standpoint of heat transfer, then, it is indicated that the skip vehicle 
is inferior to both the ballistic and glide vehicles. That is to say, 
proportionately more coolant of one form or another would be required to 
protect the skip vehicle than would be required to protect ballistic or 
glide vehicles of the same range. The skip vehicle has other disadvantages 
as well. Certainly one of the most serious of these is the very high 
lateral loads (see fig. 3) that the vehicle would be required to with- 
stand during a skip from the earth's atmosphere. These loads, coupled 
with simultaneous high thermal stresses (due to high convective rates), 
would require the structure to. be stronger and, consequently, heavier 
than that of a comparable glide vehicle. 3 For these and other reasons 
concerned with problems of stability, control, and guidance, the skip 
vehicle is thought to be the least promising of the three types of. hyper- 
velocity vehicle considered here. 

In essence, then, the preceding study has indicated that the ballistic 
vehicle exhibits the possibility of being relatively efficient for hyper- 
velocity flight by virtue of the fact that aerodynamic heating can be 
markedly reduced through the artifice of using blunt, high-pressure-drag 
re-entry shapes. The disadvantage of using the relatively inefficient 
ballistic trajectory is counterbalanced by this advantage which tends to 
keep initial mass down by reducing coolant mass. The glide vehicle 
appears promising for hypervelocity flight because it has, coupled with 
the relatively high efficiency of the glide trajectory, the possibility 
of radiating a large fraction of the heat absorbed by convection. 

Up to this point we have considered the performance efficiency of 
the several types of hypervelocity vehicle by comparison with each other. 
It is of interest now to compare, insofar as is possible, the efficiency 
of flight of these vehicles with that of lower speed, more conventional 
type aircraft. 


3 Added weight means, of course, added coolant (see, again, eq. (6l)) 
which, in turn, means added weight. The performance efficiency of vehi- 
cles is reduced accordingly - indeed one can easily demonstrate that ulti- 
mately the coolant is being added to cool coolant. This situation must 
obviously be avoided. ; 
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Comparison of Hypervelocity Vehicles 
With the Supersonic Airplane 


In the analysis of powered flight it was found that the basic per- 
formance equation for hypervelocity vehicles could be written in a form 
analogous to the Breguet range equation. Thus, according to equations ( 7 ) 
and (9), we have for both hypervelocity and lower speed vehicles that 

R - (jj|) (84) 

where it is understood that the effective quantities are the same as the 
actual quantities in the case of the lower speed, more conventional air- 
craft. Wow let us consider the product (L/D) e IeVe* Taking first the 
supersonic airplane we assume flight at a Mach number of 2 and a maximum 
lift-drag ratio of 5* The product I e V e for a ram-jet or turbojet is 
not likely to exceed about 4.4x10 s feet. 4 * The product (L/D) e I e V e is 
then 22xl0 6 feet for the airplane. Now let us compare these quantities 
with the corresponding quantities for a ballistic vehicle and let us pre- 
sume that the range will be half the circumference of the earth. In this 
event, the effective lift-drag ratio for the ballistic vehicle is 2jt (see 
fig. 2) which is slightly greater than that for the airplane, while the 
effective velocity is just half the satellite velocity, or 13,000 feet 
per second. Let us again assume that the effective specific impulse is 
225 seconds. In this case, the product of I e V e is 2.9x10 s feet which 
is about two thirds of that for the airplane. The product (L/D) e IeVe is 
about 18x10 s feet which is less than, but certainly comparable to, that 
for the supersonic airplane. Thus we have our first suggestion that the 
hypervelocity vehicle is not necessarily an inefficient type vehicle for 
long-range flight. 

In order to pursue this point further, a performance efficiency 
factor (see eq. (10)) defined as 



has been calculated for ballistic and glide vehicles for I e = 225 sec- 
onds, and ranges up to the circumference of the earth. The corresponding 
quantity E has been calculated for the supersonic airplane 


4 This observation holds essentially for any air-breathing engine - 
note that the maximum value of I e V e is simply the product of the thermal 

efficiency (taken as 0.3 - see, e.g., ref. 6) and the specific heat 

content of the fuel (taken as 1^. 6x10 s feet for gasoline-type fuels). 
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(l e V e = 4 . 4x10 s feet) for several lift-drag ratios. The results of these 
calculations are presented in figure 10 and we observe, as our example 
calculation suggested, that both the ballistic and glide vehicles compare 
favorably with the supersonic airplane for ranges in the neighborhood of 
and greater than half the circumference of the earth. The glide vehicle 
is again superior to the ballistic vehicle at lift-drag ratios in excess 
of 2 and, as a result, it compares favorably with the airplane at shorter 
ranges than the ballistic vehicle. 

It should be kept in mind, of course, that m.f may be substantially 
greater than mp, the mass of the pay load. This point takes on particular 
significance with regard to expendable vehicles where m^/mp is perhaps a 
better measure of cost than m^/mf. Thus, noting that mj_/mp=(mj_/mf ) (m^/mp) , 
and recognizing that mp/mp is probably lowest for the ballistic vehicle, 
we anticipate that the ballistic vehicle would appear to better advantage 
than shown in figure 10. 


CONCLUDING REMARKS AND SOME DESIGN CONSIDERATIONS 
FOR GLIDE VEHICLES 


During the course of this study it has been indicated that ballistic 
and glide vehicles can be operated at hypervelocities with the reasonable 
assurance that problems of aerodynamic heating can be largely alleviated 
by proper design. Skip vehicles appeared substantially less promising 
in this as well as other respects. It was further demonstrated that on 
the basis of equal ratios of initial to final mass, the long-range hyper- 
velocity vehicle compares favorably with the supersonic airplane. These 
considerations suggest that the ballistic and glide vehicles have, in 
addition to the advantages usually ascribed to great speed, the attractive 
possibility of providing relatively efficient long-range flight. 

In view of these findings, it seems appropriate as a final point to 
touch on what appear to be favorable design features of glide-type vehi- 
cles. Comparable aspects of the ballistic vehicle are not treated here 
inasmuch as they have already been considered in some detail in refer- 
ence 3 . Two categories of glide vehicle will be considered. The first 
category is made up of those vehicles whose flight through the earth’s 
atmosphere is entirely at hypersonic speeds. More specifically, a vehi- 
cle in this category is required to be stable and controllable to the 
point of high-speed impact with the surface of the earth. The second 
category includes those vehicles which are required to have acceptable 
low- speed aerodynamic characteristics (perhaps to the point of landing). 

Considering now the first category, we recognize that while the 
shape of the vehicle must be such as to provide reasonably high lift-drag 
ratios, it should also be a compact configuration designed to minimize 
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structural and, hence, propellant weight. Furthermore, it appears most 
desirable from the standpoint of aerodynamic heating to eliminate all 
surfaces that .present extremely severe heating problems. These considera- 
tions tend naturally to focus our attention on configurations free of 
the planar surfaces normally used for developing lift, and providing stable 
and controlled flight. We pursue this point, therefore, by inquiring of 
the efficiency of a body of revolution as a lifting device at hypersonic 
speeds. For this purpose, theory, including friction drag estimates, was 
employed to calculate the maximum lift-drag ratios of cones and flat 
plates (see Appendix E) and the results are presented in figure 11. It 
is seen that while the flat plate is, as would be expected, by far the 
most efficient lifting device at low supersonic speeds (l/M — >l) , it has 
relatively little advantage over the cone at hypersonic speeds (l/M— >0). 
Evidently, the body of revolution is relatively efficient for developing 
lift at hypersonic speeds. 

Now, it is clear (see eq. (78)) that just as with the ballistic 
vehicle, the nose of the glide vehicle should be rounded to alleviate 
the local heating problem. There is evidence, both theoretical and exper- 
imental, that in addition to alleviating the heating problem, rounding the 
nose may, in fact, increase the lifting efficiency of a body.- Recalling 
that for slender bodies the maximum lift-drag ratio is governed primarily 
by zero-lift drag, we recognize the validity of this statement on the 
basis of the theoretical and experimental work of reference 7 and the 
experimental work of reference 8. 

The body must, of course, be stable and controllable in flight. It 
appears feasible to provide these requirements without recourse to planar 
surfaces. Stability in pitch and yaw can be provided by a conical flare 
at the base of the body, and control in pitch and yaw can be provided by 
deflectable sections Of the surface of the body. These sections are 
located on the rear portion of the body to provide a configuration which 
is inherently stable in roll. A hypervelocity glide vehicle in the first 
category might, then, in view of these considerations, appear something 
like that shown in figure 12. 

Consider 'now a glide vehicle ' falling in category two. First of all, 
it appears most unlikely that acceptable low-speed aerodynamic character- 
istics can be obtained without using more conventional planar surfaces, 
at least to the extent of a wing. The question then is: What can be 

done to alleviate the aerodynamic heating of planar surfaces? Especially 
in this regard are we concerned with the very severe heating encountered 
at the leading edges of these surfaces. It is apparent by analogy to the 
nose of the body that the severity of aerodynamic heating at the leading 
edge of a wing can be reduced by simply rounding the leading edge. In 
fact, on comparing the theoretical results of references 9 and 10, we see 
that for the purposes of this report there is no essential difference 
between the heat-transfer rate at the stagnation point of a blunt leading 
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edge and that at the stagnation point of a comparable blunt nose as given 
by equation ( 44 ). We may anticipate, however, that rounding the leading 
edge will incur a drag penalty which will, in turn, reduce the attainable 
lift-drag ratios. This difficulty may be largely circumvented by simply 
sweeping the leading edge of the wing. The contribution to total drag 
of the drag at the leading edge is in this manner reduced approximately 
in proportion to the square of the cosine of the angle of sweep. Equally 
important, the rate of heat transfer to the leading edge may also be sub- 
stantially decreased by sweep. This possibility is suggested by the 
independence principle for cylindrical viscous flows (refs. 11 and 12) 
which applies to the components of flow normal and parallel to the axis 
of a cylinder. 5 We conclude, therefore, that in the event rounded lead- 
ing edges are used to alleviate heat transfer to wings of glide vehicles, 
or for that matter any hypersonic vehicle, these wings might well have 
highly swept leading edges. One is lead naturally to consider triangular 
or delta plan form configurations. By entirely analogous reasoning, the 
triangular plan form may also prove desirable for stabilizing and control- 
ling surfaces.. 

It is entirely possible, of course, and perhaps even desirable, 
that wing-body junctures of a glide vehicle should not be discontinuous 
but rather that the body, in effect, should be simply flattened out to 
appear more or less elliptic in cross section. In any event we see that, 
interestingly enough, the concept of sweepback may play an important role 
in reducing the heating and drag of practical configurations at hypersonic 
speeds, much as it has for drag alone at supersonic speeds. 


Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif., Dec. 10, 1954 


5 According to the independence principle and equation ( 44 ) , heat- 
transfer rates will be reduced by the cube of the cosine of the leading- 
edge sweep angle. This principle is, of course, not strictly applicable to 
viscous flow about blunt cylinders. Actually, there may be a sizable con- 
tribution to heat transfer by the component of flow parallel to the leading 
edge, and unpublished theoretical and experimental results indicate that 
the reduction in heat transfer may be more nearly proportional to the 
square of the cosine of the sweep angle. In any case, it is sufficient 
for the purposes of this report, to note that a sizable reduction in heat 
transfer is achieved by sweeping a blunt leading edge. 
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APPENDIX A 
NOTATION 

reference area for lift and drag evaluation, sq ft 
specific heat of vehicle material, ft-lh/slug °R 
drag coefficient 
lift coefficient 
skin-friction coefficient 

equivalent skin-friction coefficient, (see eq. (40)) 
specific heat of air at constant pressure ft-lh/slug °R 
specific heat of air at constant volume, ft-lh/slug °R 
drag, lh 

Naperian logarithm base 

performance efficiency factor, (see eq. (85)) 

general functional designation 

functions of AJ, (see eqs. (74) and (80)) 

acceleration due to force of gravity, ft/sec 2 

convective heat-transfer coefficient, ft-lh/ft 2 sec °R 

convective heat transferred per unit area (unless otherwise 
designated ), ft-lh/ft 2 

specific impulse, sec 

range parameter for glide vehicle (see eq. (68)) 

Stefan-Boltzmann constant for hlack hody radiation (3-7x10"" 10 
ft-lh/ft 2 sec °R 4 ) 

lift, lh 

mass, slugs 

Mach number 
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convective heat transferred (unless otherwise designated), 
ft-Ib 

distance from center of the earth, ft 
radius of curvature of flight path, ft 
radius of earth, ft 
range, ft 

distance along flight path, ft 
surface area, sq ft 
time, sec 

temperature (ambient air temperature unless otherwise spec- 
ified) , °R 

velocity, ft/sec 

ratio of velocity to satellite velocity 

velocity of satellite at earth's surface (taken as 25*930 
ft/sec) 

weight, lb 

vertical distance from surface of earth, ft 

angle of attack, radians unless otherwise specified 

constant in density-altitude relation, (22,000 ft) 1 (see 

eq. (15)) 

ratio of specific heats, Cp/C v 

semivertex angle of cones, radians unless otherwise specified 
increment 

lift-drag efficiency factor, (see eq. (C27)) 

angle of flight path to horizontal, radians unless otherwise 
specified 


P 


function of Mach number, (see eq. (E7)) 
air density, slugs/cu ft (p 0 = 0.0034) 
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nose or leading-edge radius of body or wing, ft 

partial range , radians 

total range, radians 

remaining range ($ - <p) > radians 

Subscripts 

conditions at zero angle of attack 
conditions at end of particular rocket stages 

conditions at point of maximum average heat-transfer rate 
average values 

conditions at point of maximum local heat-transfer rate 

convection 

effective values 

conditions at entrance to earth's atmosphere 
conditions at exit from earth's atmosphere 
conditions at end of powered flight 
initial conditions 
local conditions 

ballistic phases of skip vehicles 
total number of rocket stages 
pressure effects 
pay load 

recovery conditions 
radiation 

stagnation conditions 
total values 


wall conditions 
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APPENDIX B 

SIMPLIFYING ASSUMPTIONS IN THE ANALYSIS OF THE GLIDE TRAJECTORY 


The assumption of small deflection angle (0 < < l) was used through- 
out the study of the glide trajectory. In addition, equation (28) was 
simplified on the assumptions that 

(L/D)0 < < 1 (Bl) 


and 


■g do ^ L/P dV 2 
ds 2 ds 


(B2) 


The extent to which these assumptions are permissible can he checked by 
deriving an expression for (l/d) 0 and examining its variation over a 
range of trajectory parameters. 


From equations ( 30) , (33), and (3*0 the altitude of any point in a 
glide trajectory is found to he 


1 

In 

k 

1 - (1 - V f 2 )e J 

+ In 

1 - V*r*r 0 


0 ' 

(1 - V f 2 )e J __ 

V 2 r=r 0 

N— mm* 



(B3) 


By retaining the assumption of small inclination angle, whereby Q ~-dy/ds, 
and recalling that J = (2s/r 0 )/(L/D) , we find the inclination angle by 
differentiating equation (B3) . Performing this operation and making use 
of equation (30) reduces the expression for (L/D )0 to 


2 

1 

2.105X10 

0r o 

-1 - (l-V f 2 )e J - 

' V 2 


-3 


(b4) 


Since V becomes "very small toward the end of the trajectory, it is 
apparent from equation (B4) that the assumption of small (L/D )(9 cannot 
be justified. in this portion of flight. The problem then is to determine 
the conditions under which (L/D)0 remains negligibly small over the 
major part of the trajectory. 
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With the aid of equations (30) and ( 36 ), equation (B4) can he modi- 
fied to the following form 


l = 1 - 42 ln ( 
$ $ 


1 


1 - 


g/Pr 0 . 
,(L/D)eJ 


(B5) 


For given values of l/d and total range 0 , equation (B5) determines 
the fractional part of the total range which corresponds to a given value 
of (L/D)0. Since the deflection angle is always increasing, we can there- 
fore determine the portion of the total range through which (l/d )0 remains 
equal to or less than a given value. A computation of this nature was 
performed for a value of (l/d )e < 0 . 05 , and the results are presented in 
figure 13. From this figure we can see that except for short ranges 
and large lift-drag ratios, (l/d )0 (as well as 0 ) remains at a value less 
than 0.05 for better than 90 percent of the total range. 

The second assumption, equation (B2) , can also be verified from the 
results of the analysis. By differentiation of equation (B4) we find that 


V 2 ^ = 


4v c 


(1 - V f 2 )e' 


ds pr 0 2 (l/D) : 


V 2 


(B 6 ) 


while differentiation of equation ( 30 ) yields 


I cw»> fr = " 


Vs 2 (l - V f 2 )e J 


(B7) 


Dividing equation (B 6 ) by equation (B7) , and making use of equation (B4) , 
we find that 


V 2 


d0 

ds 


\ (l/d) 


dV 2 

ds 


— — = 2 — t(L/b)0] 

(L/D) (L/D) 2 


(B 8 ) 


By comparing equation (B 8 ) with the previous results obtained for (l/d) 0, 
(fig. 13) we can readily see that this assumption, equation (B2) , is 
actually less stringent than the previous one for practical values of 

L/D. 



NACA TN 4o46 


39 


APPENDIX C 


THE RELATION BETWEEN 


Cjp’S 

— — AND 

c d a 



FOR CONICAL MISSILES 


The lift and drag coefficients for slender cones at small angles of 
attack can he expressed in the following manner: 


C L = 2a 
C D “ C D 0 + 

By dividing the first of the preceding 
obtain an expression for the lift-drag 

L _ = ^L _ ^L 

D C D C Dq + aC L C Dq + C L 2 /2 

It can be shown that equation (C3) has a maximum value when 

C Dq = C L 2 /2 (C4) 


(Cl) 

aC L (C2) 

equations by the second, one can 
ratio 


whereby 


(CD) (L/D) r - ’ 


'max 


(C5) 


From equations (Cl) through (C5), the maximum value of the lift-drag 
ratio can be expressed in the following ways: 

= JL = _L = _JL_ (c 6) 

VVmax 2a Cl J2Cj) q 

The drag coefficient at zero angle of attack appearing in equa- 
tion (C2) can be broken down into its component parts to yield 



(C7) 



4o 
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where (Cn ) is the zero-lift pressure drag coefficient and C F is the 

u O p u 

zero-lift skin-friction coefficient based on wetted area. The skin- 
friction coefficient, Cf q , in equation (C7) can be related to the equi- 
valent skin-friction, Cp‘, (see eq. (1+3)) by considering average conditions 
over the. surf ace of the cone. Equating the friction drags as determined 
from free-stream and local average conditions, it is found that 



/ \ (Pz)av( v l) a y 2 

\ F Vav PV 2 


(C8) 


By referring to local average conditions on the body surface, the expres- 
sion for Cp' , equation (1+0), can be written as 


Cp« 



(Pz ) av ) ay 
“pV 


(C9) 


Comparing equations (C8) and (C9) it is apparent that 

, (Vi), 


: F - c F 


9 av 


V 


(CIO) 


For slender shapes at hypersonic speeds, the local velocity does not 
differ appreciably from the free-stream value. Also, for small angles 
of attack, the skin-friction coefficient should remain fairly constant. 
Consequently, equation (CIO) can be written as 


Ctt = Ctp = Ctp t = constant 


(Cll) 


and equation (C7) then becomes 


Cp‘S 

C D 0 " ( C D 0 )p + A 


From equations (C5) and (C12) it can then be shown that 


/Cj 

\c D A 


(l/d) 


1 

2 


max 


1 - 


(C Do )i 


r 


(C12) 


(C13) 
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From the Newtonian impact theory, the zero-lift pressure drag coef 
ficient for slender cones at hypersonic speeds can he expressed as 


(CD 0 ) p = 

where 6 is the semi vertex angle of the cone, 
for slender cones 

S I 
A " 6 


(Cl4) 

By further noting that 

(C15) 


equation (C12) then becomes 


^o 


25 



(C16) 


For a given value of Cp* it can be shown that equation (Cl6) has 
minimum value when 1 


6 = 


^opt - 



1/3 


(C17) 


whereby, at 6 = 6 Qpt 

* 6S opt 2 = 3(C„ 0 ) p (C 18 ) 


Obviously, then, the highest value of maximum lift-drag ratio (eq. ( C6) ) 
will be attained by the cone with the semivertex angle given by equa- 
tion (C17) . By substitution from equation (Cl8) into equation (C 6), the 
opt imum value of maximum lift-drag ratio is 



l ( 4 V /a 

2^ Vc F V 


1 

2*^ & 0 pt 


(C19) 


By further substituting the expression for minimum zero-lift drag 
coefficient, equation (Cl8), into equation (C13), the following relation, 
corresponding to the condition of optimum maximum lift-drag ratio, is 
obtained: 


27 
>a A 


C_F 

c D A y(L/D) 


max 


1 

3 


(C20) 


1 The remaining analysis assumes Cp* constant. Although this is 
definitely not the situation in real flight, this analysis provides an 

"order of magnitude" estimate of pertinent parameters. 
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With the aid of equations (Cl4) and (Cl6), equation (C13) can also 
be expressed in the following form, corresponding to any maximum lift- 
drag ratio including the optimum value: 


/ ys\ 

v^Ml/dW 



(C21) 


From equations (C20) and (C21) it can readily be seen that in the case of 
the optimum (L/D) max 


Cp' 

2$opt 3 


= 2 


(C22) 


from which it follows directly that 



(C23) 


so that equation (C2l) may be written as 


'C F 'S\ 

' C » A A L /D) max 



(C2^) 


With the aid of equations (Cl6) and (C23)> the expression for any 
(L/D) max (eq, (C6)) can be shown to be 


D/max 


1 

28 

1 + 2| 


3 n 

u 

1/2 


(C25) 


and it follows directly from equation (C19) that the ratio of (L/D) max 
for any cone to that for the optimum cone is 
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11 



(C2 6) 


where T} is defined as the "lift-drag efficiency factor." By substitu- 
tion from equation (C24), the ratio of (L/D) m ax to the optimum value can 
then be expressed in terms of (Cp'S/CpA) as follows: 





c F «sy /6 

c D Ay 


(C27) 


The dependence of r\ on Cp f S/CpA is shown in figure It should 

be noted, however, that for small values of q the assumption of slender 
cones will be violated, although the results as shown will be qualitatively 
correct in that C^S/C^A will become exceedingly small for low values 
of (L/D) ma x regardless of body shape. 
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APPENDIX D 


COMPUTATION OF HEATING ASSOCIATED WITH ROCKET VEHICLES 
Radiation of Heat From Glide Vehicles 


From equation (72), the maximum time rate of total heat input to 
the glide vehicle can be expressed as 

/dQ\ _ av\ = m § v S fSll£] (Dl) 

Wmax V dt /max l/d) \ C V a J 


The rate of heat radiation from the vehicle can be expressed by the fol- 
lowing standard relation 


s TT = Kr4s 


(D2) 


Using equations (Dl) and (D2) , the requirement for continuous radia- 
tion of all convective heat input to a surface at a temperature 2000 R 
(1540° F) can be expressed as 


< i.2o (D3) 

s(l/d) c d a 


If a value of (L/D) max = 6 is assumed, values of the parameter C F *S/C D A 

opt 

and cone angle, 5, can be determined as a function of actual (L/D) max 
from the analysis given in Appendix C. A vehicle weight of 5000 pounds 
with a maximum diameter of 3 feet is assumed whereby equation (D3) can 
be evaluated for various (l/d) giving the results in the following table: 


l/d 

CF’S/CpA 

6, 

deg 

s, 

sq ft 

mg Cp'S 

S(l/D) CpA 

6 

0.333 

2.75 

147 

1.89 

4 

.0600 

6.73 

6 0.2 

1.24 

2 

'.00710 

14.3 

28.6 

.619 

1 

.000890 

29.6 

14.3 

.312 

1/2 

.000115 

55.8 

8.55 

.135 
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We see, therefore, that at surface temperatures of 2000° R and for an 
l/d of 4, this glide-type vehicle can radiate heat at a rate equal to or 
greater than the maximum convective heat rate. 


Radiative and Convective Heat Transfer Associated 
With Skip Vehicle 


In this section the problem is to determine the extent to which 
heat absorbed by a skip vehicle in the first skipping phase, can be 
reradiated during the subsequent ballistic phase. The quantity of heat 
absorbed in the first skipping phase has - already been obtained in the 
heating analysis, (eq. (59) for n = l) 


Qi _ 1 ( 

i raVf 2 " 5 C » A \ 



(DU) 


where the total heat absorbed throughout the entire trajectory is 


Or i 

| mv f 2 = 2 "CdA 


(D5) 


In order to determine the heat radiated, three quantities must be 
determined: 

1. Temperature of the vehicle at the start of the second ballistic 

phase 

2. Temperature of the vehicle at the end of the second ballistic 

phase 

3« The time duration of the second ballistic phase 

To determine the first of the above quantities, we employ the rela- 
tion for heat absorbed 


Qi = cWgAT (D6) 

where c is the specific heat of the material, W e is the effective 
weight of material absorbing heat, and AT is the temperature rise during 
the first skip. If it is assumed that 1/3 of the missile weight will 
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absorb heat, equation (D6) becomes 

AT - i (t) < dt> 

where m is the total mass of the vehicle. It is assumed that the 
material has a specific heat of 0.11 Btu/lb °R. If it is also assumed 
that the temperature at the start of the first skip is 500° equa- 

tion (D7) becomes 

T eX2 = 500 + l.lXlO" 3 ^ (D8) 

which defines the temperature at the beginning of the second ballistic 
phase. 

To find the temperature at the end of the second ballistic phase, 
we equate the radiant heat-transfer rate from the body to the rate of 
heat loss in terms of the temperature drop of the body 

-kT 4 S dt = cW e dT (D9) 


This expression can be integrated between limits from beginning to end 
of the second ballistic phase to yield 


T 


en 2 


( 7. 95X10" 15 ) St + 


u ex 2 


(D10) 


for a vehicle weight of 5000 pounds (effective absorbing weight of 1667 
pounds) where T en2 as the temperature at the end of the second ballistic 
phase and t is the total time of the second ballistic phase. The total 
heat lost by radiation can now be expressed in terms of the temperature 
drop as 


Qr 2 - (^ex 2 ~ ^en 2 )^e c 
or 


Qr 2 - 1. ^lxl0 5 (T e x 2 “ T e n 2 ) 


(Dll) 
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The time of flight in any ballistic trajectory can be shown to be 


t = 


2r r 


1 - cos 


I) 


3/2 


(1 - i 2 ) 


tan Qf + 


J 1 - i 2 


tan 


-1 


(TrH tan t) 


J (D12) 


vhere 


tan Qf 


1 - 


sin | + tan 0fcos ^ 


— = tan x | 
2 ' 


sin 6 fcos Qf\ 

=k - COS 2 0f 


The foregoing relations were applied to a computation of the radia- 
tive cooling of a missile weighing 5000 pounds and traversing a total 
range of 3440 nautical miles (0 = 1.0). Values of Qf were obtained in 
the motion analysis , whereas values of Cp'S/c^A and S obtained in the 
previous calculation with regard to the glide missile will apply to this 
case also. The computations are summarized in the following table. Note 
that the case of L/d = l/2 is essentially the ballistic vehicle (see 

fig. 9). 


l/d 

Of , 

deg 

Vf 2 

QlXlO -5 

Qi/Qt 

sec 

^ex2-’ 

°R 

Ten 2 > 

°R 

% 2 xl0" 5 

Qr 2 /Qi 

6 

12.5 

0.275 

3115 

0.135 

213 

2710 

1490 

1725 

0.554 

4 

17.0 

.315 

1470 

.258 

335 

1542 

1323 

316 

.211 

2 

24.0 

• 525 

549 

.575 

395 

889 

885 

27 

.049 

1 

27.5 

.620 

122 

.853 

247 

587 

585 

3 

.023 

1/2 

30.0 

.650 

19 

.985 

80.4 

514 

514 

0 

0 


We see, therefore, that the quantity of heat which must be absorbed 
by this skip vehicle decreases rapidly with decreasing lift-drag ratio. 

The quantity of heat which must be absorbed by a ballistic vehicle 
(L/D = l/2) is almost negligible compared with the quantities associated 
with vehicles with an I/D = 2 or greater. Comparison of the heat absorbed 
in the first skipping phase with the heat radiated in the second ballistic 
phase indicates no appreciable advantage is obtained due to radiation for 
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values of L/D = 2 and lower. To be sure, this situation could be sub- 
stantially altered (near L/D = 2) by allowing the surface temperatures 
to reach higher values during the skip; however, it seems most unlikely 
that the net heat absorbed by the skip vehicle could ever be reduced to 
the low value of the ballistic vehicle for any reasonable surface tern- 
perature. 
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APPENDIX E 


DETERMINATION OF THE RATIO OF (L/D) max FOR 
A FLAT PLATE TO (L/D) max FOR A CONE 

In Appendix C the expression for the optimum (L/D) max for a. slender 
cone was developed with the assumption of large Mach numbers and a constant 
value of C F '. A simple analysis will now be presented whereby It is 
shown that equation (C19) will apply, for the most part, throughout the 
Mach number range and for a variable Cp* . Inasmuch as the lift coef- 
ficient (eq. (Cl)) is essentially independent of Mach number and cone 
angle, modification of the results obtained in Appendix C (eq. (C19)) 
will occur only through the evaluation of the zero-lift drag coefficient 
(eq. (C6)). 

The variation with Mach number of the zero-lift pressure-drag coef- 
ficient, (Cp 0 )p, can be represented as 

(Cn ) = 2f (M) sin 2 S (El) 

o p 

while Cp* is assumed to vary with cone angle in the following manner 

C F « = C F » 5 m (E2) 

r r m=o 

where m is chosen to give the required variation of Cf' with &. By 
an analysis paralleling that presented in Appendix C, it can be shown 
that the optimum value of the lift-drag ratio assumes the following form 



Now the values of f(M) should not differ greatly from 1 while values of 
m should not differ greatly from 0. It follows then that values of 
(L/D) obtained from equation (E3) will not differ appreciably from those 
those obtained from equation (C19) . Furthermore, the coefficient C F * 
will not differ appreciably from C F for the slender shapes and small 
angles of attack under consideration. Equation (E 3 ) then reduces to 
the form 
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/l\ = _i_ /V\ 

^ ' max 2 n/3 'PfJ 


(E4) 


which corresponds to equation (C19) and is also a suitable approximation 
for all supersonic Mach numbers. 

For a flat plate in supersonic flow, the lift and drag coefficients 
can be expressed as 


C L = -=== + 2a 2 (E5) 

s/m 2 - 1 

C D = C D 0 + ( E 6) 


- which follows easily from the equation for pressure coefficient given by 
Ivey and Cline (ref. 13) . By noting that 


C Dq = 2C I 


the lift-drag ratio can be written from equations (E5) and (E6) as 


2Aa + a 2 


D 


Ctt + 2Ao o.. + ex'" 


(E7) 


where 


A = 




Equation (E7) is found to have the following maximum value 


/max a + 


aA 


2aA 


(E8) 


2(A + a)(2A + a) (2A + a) 


which occurs at values of the independent variables given by the maximum 
condition 
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Cjr = a 2 


2A + 


a" 


2(A + a) 


(E9) 


From equations (e4) , (E8) , and (E9) , and recalling that A = (M - l) x ^ 2 
the ratio of (L/D) max for a flat plate to (L/D) max for a cone can be 
determined as a function of Mach number and Cp. The variation of this 
ratio with Mach number for given values of Cp 1 has been presented in 
figure 11, where it was found convenient to use the inverse Mach number as 
the abscissa in order to illustrate conveniently the behavior at hyper- 
sonic speeds. 
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Figure /.— Variation of velocity with incidence angle for various values 

of range of ballistic vehicle. 
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Range parameter, & 


i 


Figure 2- Variation of effective lift -drag ratio with range 
for optimum ballistic vehicle. 
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Figure 3- Trajectory of the first skipping phase for a skip vehicle with 
a lift- drag ratio of 2 and a total range of 3440 nautical miles 

<<p -- // 
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Figure 5- Variation of effective lift- drag ratio with range for various 
values of aerodynamic lift- drag ratio of skip vehicle. 
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Range parameter, <P 


Figure 7- Variation of effective lift- drag ratio with range for various 
values of aerodynamic lift-drag ratio of glide vehicle. 
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Figure 8.- Variations of f^(AJ) and F b (AJ) with AJ. 
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Figure 9- Variation of mass ratio with range for various lift -drag ratios of hyper velocity vehicles. 
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Figure 13.- Portion of range where (L/D) 0 - 0.05 as a function of range for various 

values of lift -drag ratio. 
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